Remarks On The Spherical Scalar Field Halo In Galaxies 
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£N| ' Abstract 



o 



Matos, Guzman and Nunez proposed a model for the galactic halo within 



q-i the framework of scalar field theory. We argue that an analysis involving the 

full metric can reveal the true physical nature of the halo only when a certain 
condition is maintained. We fix that condition and also calculate its impact on 
observable parameters of the model. 



> 

f^) ■ One of the outstanding problems in modern astrophysics is the problem of 

dark matter which is invoked as an explanation for the observed flat rotation 
■ curves in the galactic halo. Doppler emissions from stable circular orbits of 

' neutral hydrogen clouds in the halo allow the measurement of tangential ve- 

locity Vt g (r) of the clouds treated as probe particles. According to Newton's 



laws, centrifugal acceleration v\ jr should balance the gravitational attraction 



GM(r)/r 2 , which immediately gives v\ = GM(r)/r. That is, one would expect 



"tgl 

res'v tg 

km \ a fall-off of Vf g (r) with r. However, observations indicate that this is not the 

case: v tg approximately levels off with r in the halo region. The only way to 
interpret this result of observation is to accept that the mass M(r) increases 
linearly with distance r. Luminous mass distribution in the galaxy does not 
follow this behavior. Hence the hypothesis that there must be huge amounts of 
nonluminous matter hidden in the halo. This unseen matter is given a technical 
name dark matter. 

Despite the fact that the exact nature of dark matter is as yet unknown, 
several analytic halo models exist in the literature including those provided 
by scalar-tensor theories (see for instance [1]). In particular, the scalar field 
model first proposed by Matos, Guzman and Nunez [2] has received considerable 
attention. It is important to note that the authors primarily constructed an 
exact solution of Einstein's field equations sourced by a scalar field that provides 
a density profile of 1/r 2 together with other appealing features of the metric 
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functions. As a particular application, they sketched a plausible interpretation 
of the halo dark matter problem. The problem being important in itself, we 
think that the interesting relativistic central feature of the solution, namely, a 
non-Newtonian halo, must be well grounded. The purpose of the present Brief 
Report is to fix the condition under which it is possible. In addition, we work 
out its impact on observable parameters. 

It is to be mentioned that the solution in [2] has been criticized because of 
its singular behavior at the origin [3] , but this singularity is not peculiar to that 
solution alone; there are other viable halo models in the literature that also 
possess such singularity (see for instance [4] ) . While we are here interested only 
in the outer reaches of the galaxy (very far from the origin) , the authors of Rcf . [2] 
and the associated research group have addressed this inner region singularity 
too. They obtained several new results under the scalar field dark matter model 
in galaxies: Solution with axial symmetry including the inner region [5], time- 
dependent spacetimes [6] , the full non- linear Newtonian evolution after the turn- 
around point [7], time evolution of density fluctuation [8], collision properties 
of two structures [9] and so on. While they obtained constraints arising out of 
cosmological scale considerations, we believe that it is also useful to ascertain 
the constraint arising out of the local scale, which would clarify the relativistic 
nature of the spherically symmetric model of the halo under consideration. 

Using the flat rotation curve condition [10], Matos, Guzman and Nunez 
obtain the spherically symmetric static solution for the galactic halo as follows 
(G = c = 1, unless specifically restored): 

ds 2 = -B(r)dt 2 + A(r)dr 2 + r 2 (d9 2 + sin 2 6d(f> 2 ) (1) 

B(r) = B r l 
A -I 2 

A{r) = A + D{A-l 2 )r-( l + 2 ) (2) 




0M = a/ — mr + o (3) 



V V = (4) 

where D is an arbitrary constant of integration, <j> and V are the scalar field and 
potential respectively. The parameter I = 2(v tg /c) 2 , B > is another constant. 
Observations of the frequency shifts in the HI radiation show that, in the halo 
region, v tg /c is nearly constant at a value 7 x 10~ 4 [11]. Thus, in what follows, 
we take I ~ 1CT 6 . 

Note that we can rewrite A(r) in the standard Schwarzschild form 



A(r) 



^ 2m(r) 



(5) 



which is often convenient and will be useful later while discussing the observa- 
tional parameters. Such a form has the advantage that it immediately reveals 
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not only the mass parameter m(r) but also shows that the proper radial length 
is larger than the Euclidean length because r > 2m(r). This inequality, which 
is essential for signature protection, dictates that A(r) > 1. This is a crucial 
condition to be satisfied by any valid metric. 

Now, for the sake of simplicity, Matos, Guzman and Nunez choose D = 0, 
but this is not the best choice because it makes the metric component A < 1. As 
a consequence, whatever results follow from the reduced metric should be taken 
with caution. For instance, the stresses exhibit a density profile p < 0, meaning 
violation of Weak Energy Condition (WEC) and furthermore lead to u < — 1 
(see below), meaning repulsive gravity in the halo, contradicting observational 
facts. But these are actually not the true features of their model. To see the 
true picture, it is necessary to calculate the relevant quantities with D ^ 0. 

We find the density and pressure profiles in the rest frame of the fluid as 

1 r-( 4 +') [D(l 3 + I 2 - Al - 4) + l 2 r 2+l ] 

(6) 



P 8tt I 2 - A 

1 r-( 4 +') [D(l 3 + l 2 -Al-A)- l(A + l)r 2+l ] 
Pr ~ 8n I 2 - 4 (7) 



1 r-( 6 +') [D(l 3 + I 2 - Al - 4) + l 2 r 2+l ] [{r 2 - 1)1 - 2(r 2 + 1)] 
Pt ~^ 4(^4) (8) 

where p is the energy density, p r is the radial pressure and p t are the transverse 
pressures. 

Matos, Guzman and Nunez conclude that their model has huge pressure over 
density and thus it is non-Newtonian. We wish to emphasize that the role of non- 
zero value of D is crucial not only for avoiding repulsive gravity (as alluded to 
above) but also for arriving at a correct conclusion about the relative strengths 
between pressure and density. For instance, let us take D = 1. In the distant 
halo region, we can take, typically, r ~ 100 — 300 Kpc and with I ~ 10~ 6 , we find 
the numerical values to be p ~ 10~ 9 and p r ~ 10~ 9 , which means that they are 
of the same order. But on the other hand, p r + 2p t ~ 10~ n => p r + 2p t ~ 10~ 2 p, 
which indicates that total pressure is roughly one hundred times less than the 
density. However, if we take D = 10~ 5 , we find that p r + 2p t ~ 10 3 p. If we 
keep on decreasing the value of D further (but never exactly to zero for reasons 
stated above), we see that the total pressure dominates more and more over 
density reinforcing the non-Newtonian nature. 

The next question is: How far can we go on decreasing Dl We notice the 
following interesting scenario. When D = 10~ 7 , we find p r + 2p t = 9 x 10 5 p, 
which leads to w = Pr tp Pt = 3 x 10 5 (attractive gravity) as shown in Fig.l. 
This is the extreme possible non-Newtonian halo in the scalar field model under 
consideration. The reason is this. If D = 10~ 8 , we find that u> > up to 
r = r = 200 Kpc (attractive gravity) and becomes u < — 1 after r = r 
(repulsive gravity). At r = r , there is a singularity in u. This value of D 
represents a transition from attraction to repulsion as shown in Fig. 2. When 
D < 10~ 9 , we find that p < 0, ui < — 1 (repulsive gravity), which share the 
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woes that follow also from the choice D = (Figs. 3 and 4). These show that 
we can not decrease D below 10~ 7 , that is, we must have D > 10~ 7 . This is 
the condition that must be maintained in order to have a non-Newtonian halo. 

The pressures are anisotropic, as is evident from Eqs.(7,8), which is a good 
feature of the solution from the point of view of exterior matching. Note that 
the solution can not be matched to the Schwarzschild exterior metric at the 
boundary of the halo if the pressures were isotropic [12]. It can be further 
verified that p > 0, p + p r > 0, p + p r + 2p t > for D > 10~ 7 ; so we can say 
that the halo matter is not exotic because the standard energy conditions are 
satisfied everywhere. Therefore, we expect attractive gravity in the halo. To 
confirm it, we follow the prescription by Lyndcn-Bcll, Katz and Bicak [13], and 
find that the total gravitational energy is indeed negative: 



rr 2 i 

Eq = 47r / [1 - Az]pr 2 dr < 0, 

Jri 



(9) 



due to the fact that p > 0, 1 — A? < and r2 > r\. This prescription has been 
very useful in the case of scalar field wormholes too [14-16]. 

Certainly, the scalar field model corresponding to D = 1CP 7 is highly non- 
Newtonian because p r + 2p t ~ 10 6 p. As a result, a purely Newtonian definition 
of mass, viz., M{r) = 4ir f pr 2 dr does not apply. However, incorporating the 
pressure contribution, the dynamical mass in the first post Newtonian order 
becomes 

M pN (r) =4tt J(p + p r + 2p t )r 2 dr = 10 6 M(r), (10) 

which clearly reflects the non-Newtonian nature of the model in terms of masses. 

We next focus on the observable parameters expected in this non-Newtonian 
halo. Whatever be the analytic model for it, there must be a way to contrast 
its predictions with actual measurements. The key point is that one does not 
directly measure the metric functions but indirectly measures gravitational po- 
tentials and masses from rotation curve and lensing observations. Faber and 
Visser [17] have shown how, in the first post Newtonian approximation, the 
combined measurements of rotation curves and gravitational lensing allow in- 
ferences about the mass and pressure profile of the galactic halo as well as its 
equation of state. 

The usual techniques for obtaining the potential for rotation curve (RC) 
measurements yield a pseudo-potential (See Ref.[17] for details): 

$rc = $ + &n, (11) 
where $ at is the Newtonian potential, $ = \ In B and a pseudo-mass 

m RC = r 2 $'(r) w 4tt J ' (p + p r + 2p t )r 2 dr. (12) 

Faber and Visser also define the lensing pseudo-potential as 

<&(r) 1 f m(r) 



$lons = 
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and a pseudo-mass mi ons obtained from lensing measurements as 

mi ons - \r 2 ^'{r) + ^m(r). (14) 
The first order approximations of Einstein's equations yield 

w 4^2[2<nsW - m' RC (r)} (15) 

47rr 2 (p r + 2pt) » 2[m' RC (r) - m[ ens (r)\ (16) 

where the right hand sides denote pseudo-density and pseudo-pressures. Fur- 
thermore, Faber and Visser define a dimcnsionless quantity 

w{r) = Pr + 2p t w 2mkc-mU , (17) 
U 3p 3 2m( ens -r^ c 1 J 

The pseudo quantities on the right hand side of Eqs.(ll)-(17) are actual 
observables from the combined measurement. If the observed pseudo profiles 
reasonably match with the analytic pseudo profiles coming from a priori given 
metric functions, one can say that the solution is physically substantiated. Oth- 
erwise, it has to be ruled out as non- viable. The impact of a small non-zero D on 
the analytic pseudo profiles can now be computed. For the extreme (D = 10~ 7 ) 
Matos, Guzman and Nunez solution, these work out to leading order in r as 

™ncM = | » KT 6 r (18) 
m lo nsM « « (19) 

2(m' RC - m[ cns ) » Kl *~ l _~* ] « 10- 6 . (20) 
The dimensionless parameter u to all orders in r with no restriction on D is 

1 j ~ 3 2m; cns - m' KC 3[D(I 3 + Z 2 - 4Z - 4) + Z 2 r 2 +«] ' [ ' 

which yields w«3x 10 5 for D = 10~ 7 within our chosen range, r ~ 100 — 300 
Zfpc. Note that if we straightaway put D = in Eq.(21), we get ui(r) < — 1, 
conveying a completely wrong physical conclusion. 

The pivotal result of the present article is that D > 10~ 7 and not D = 0, as 
discussed above. Of course, the lowest limit on D is small and it is quite tempting 
to set it exactly to zero. But the price for it is that one gets a completely 
wrong picture of the halo. We have analyzed the model taking into account 
only the lowest value of D. Similar analysis can be carried out with other 
values of D as well respecting the suggested lower limit. We can say that, by 
and large, the conclusion of Matos, Guzman and Nunez about non-Newtonian 
nature of the halo is right provided the restriction on D is maintained. With 
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this restriction in place, their model can indeed be a physically viable one. If 
combined measurements follow the pattern as indicated in Eqs.(18)-(21), we 
would say that the model is observationally supported. However, given the 
present uncertainties in observation, it is yet too premature to say so. 

We are deeply indebted to Guzel N. Kutdusova for her assistance at SSPA 
and BSPU where the work was carried out. 

Figure captions: 

Fig.l. Plot of u(r) vs r in which w is computed from either Eqs.(6)-(8) or 
(21) with I = 1CP 6 and D = 10~ 7 . The distance r in galactic halo region is 
taken in the range 100 — 300 Kpcs. The non-Newtonian values of ui are evident. 

Fig. 2. Plot of uo(r) vs r in which w is computed from either Eqs.(6)-(8) or 
(21) with / = 10~ 6 and D = 10~ 8 . The distance r in galactic halo region is 
taken in the range 100 — 300 Kpcs. The figure displays the transition behavior 
of oj as discussed in the text. 

Fig. 3. Plot of ui(r) vs r in which co is computed from either Eqs.(6)-(8) or 
(21) with I = 10~ 6 and D = 10~ 9 . The distance r in galactic halo region is 
taken in the range 100 — 300 Kpcs. The values of u> are negative indicating 
repulsion. 

Fig. 4. Plot of Lu(r) vs r in which u) is computed from either Eqs.(6)-(8) or 
(21) with I = 10~ 6 and D = 0. The distance r in galactic halo region is taken 
in the range 100 — 300 Kpcs. The values of u are negative indicating repulsion, 
similar to that in Fig. 3. 

References 

[I] S. Fay, Astron. Astrophys. 413, 799 (2004) 

[2] T. Matos, F.S. Guzman and D. Nunez, Phys. Rev. D 62, 061301 (2000) 
[3] U. Nucamendi, M. Salgado and D. Sudarsky, Phys. Rev. D 63, 125016 
(2001) 

[4] F. Rahaman, M. Kalam, A. DeBencdictis, A. A. Usmani and Saibal Ray, 
Mon. Not. R. Astron. Soc. 389, 27 (2008) 

[5] T. Matos, F.S. Guzman, Class. Quantum Grav. 18, 5055 (2001) 
[6] M. Alcubierre, F.S. Guzman, T. Matos, D. Nunez, L.A. Urena-Lopez and 
P. Wicdcrhold, Class. Quantum Grav. 19, 5017 (2002) 

[7] F.S. Guzman and L.A. Urena-Lopez, Phys. Rev. D 68, 024023 (2003) 
[8] F.S. Guzman and L.A. Urena-Lopez, Astrophys. J. 645, 814 (2006) 
[9] A. Bernal and F.S. Guzman, Phys. Rev. D 74, 103002 (2006) 
[10] S. Chandrasckhar, Mathematical Theory of Black Holes (Oxford Uni- 
versity Press, 1983) 

[II] A. Boriello and P. Salucci, Mon. Not. R. Astron. Soc. 323, 285 (2001) 
[12] S. Bharadwaj and S. Kar S, Phys. Rev. D 68, 023516 (2003) 

[13] D. Lyndcn-Bell, J. Katz and J. Bicak, Phys. Rev. D 75, 024040 (2007) 
[14] K.K. Nandi, Y.Z. Zhang, R.G. Cai and A. Panchenko, Phys. Rev. D 

79, 024011 (2009) 

[15] K.K. Nandi, A. Islam and J. Evans, Phys. Rev. D 55, 2497 (1997) 
[16] K.K. Nandi, B. Bhattacharjee, S.M.K. Alam and J. Evans, Phys. Rev. 

D 57, 823 (1998) 

[17] T. Fabcr and M. Visser, Mon. Not. R. Astron. Soc. 372, 136 (2006) 



6 



